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Abstract 

We discuss the properties of eigenphases of S'-matrices in random mod- 
els simulating classically chaotic scattering. The energy dependence of the 
eigenphases is investigated and the corresponding velocity and curvature dis- 
tributions are obtained both theoretically and numerically. A simple formula 
describing the velocity distribution (and hence the distribution of the Wigner 
time delay) is derived, which is capable to explain the algebraic tail of the 
time delay distribution observed recently in microwave experiments. A de- 
pendence of the eigenphases on other external parameters is also discussed. 
We show that in the semiclassical limit (large number of channels) the curva- 
ture distribution of S-matrix eigenphases is the same as that corresponding 
to the curvature distribution of the underlying Hamiltonian and is given by 

the generalized Cauchy distribution. 
72.20.Dp,05.45+b,72.10.Bg. 
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I. INTRODUCTION 



Quantum chaotic scattering has been discussed for a number of years It may occur 
in a variety of different physical situations from atoms and nuclei to disordered mesoscopic 
devices or microwave cavities. The schematic model of the system is presented in Fig.|I|. The 
cavity (internal region) is coupled to the outside world by leads. The characterictics of the 
internal motion manifest themselves obviously in the properties of the ^-matrix. This is 
seen directly using the Hamiltonian approach to the scattering B. 

Consider a simple Hamiltonian with the Hilbert space spanned by N discrete states \k > 
and M continua \c,E >: 

H = ZLi E k \kxk\+ Y.cU I dEE\c, E><c,E\ + 

(1.1) 

+9 EcU Eti I dE (W kc (E)\k ><c,E\ + h.c.) . 

Note that no continuum- continuum coupling is permitted in the model. The bound- 
continuum coupling is characterized by coupling constant g and the energy dependent matrix 
Wkc{E), where it is assumed that the columns of W are normalized to unity. If W depends 
only weakly on the energy, the corresponding unitary M x M matrix S may be expressed 
as §! 

S cd {E)=5 cd -2g 2 tY^W ck {E)(— 1— ) Wu(E), (1.2) 

ki \h, — n) n 

where TC is the effective Hamiltonian describing the motion within the bound subspace after 
eliminating the continua in a Markov approximation. In an arbitrary basis spanning the 
bound subspace it takes the form 

At 

H kl = H kl -ig 2 Y.W kc W d (1.3) 

c=l 

The question may be posed whether, for generic systems, there is a unique relation be- 
tween S-matrix properties and the type of motion inside the cavity. One way to address this 
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issue is via the semiclassical theory p] which seems, however, to be limited to a large number 
of channels, M. On the other hand, in recent experiments on scattering in microstructures 
0, M can be of the order of unity. 

A second possible way is a stochastic approach in which the Hamiltonian H and the 
coupling matrix W |2| are modeled by random matrices. A subsequent averaging over 
different realizations of H (typically for fixed W) yields statistical predictions concerning 
fluctuations of physical quantities of interest. One hopes, extending the conjecture which 
has been quite useful for bounded chaotic systems Q] that the properties of fluctuations 
are universal. For bounded, autonomous, classically chaotic systems, depending on their 
symmetries, the statistical spectral fluctuations are, generically, well represented by the 
corresponding quantities obtained from Gaussian Orthogonal (GOE), Unitary (GUE), or 
Symplectic (GSE) Ensembles of random matrices ||. 

Assuming that H matrices are drawn from one of these ensembles and making similar 
assumptions on the coupling matrix W one is forced to ask what are the properties of the en- 
semble of S-matrices. A partial answer has been obtained by Lewenkopf and Weidenmuller 
0: if if belongs to GOE and the channels are equivalent (see below) S matrix belongs 
to COE provided the coupling constant g equals to unity. This result has been recently 
generalized to all three universality classes by Brouwer || for arbitrary g. Brouwer found 
that the S'-matrices, Eq. (|1.2| ), conform with Generalized Orthogonal (GCOE) and Unitary 
(GCUE) Circular Ensembles [plj which reduce to circular ensembles of Dyson (COE, CUE), 
|TT| when the coupling to the continuum becomes ideal. 

Brouwer's result provides a direct link with another popular, random matrix theory 
based approach in which one draws directly S'-matrices, or transmission matrices, from an 
appropriate random matrix ensemble [|TQ] , |T2| , [T3|1 . The Hamiltonian approach is, in a sense, 
more general since it allows for calculation of time delays and energy-correlation averages, 
while the ensemble of S-matrices is energy independent. On the other hand, one is frequently 
interested in single-point with respect to energy statistical measures. Those may be directly 
accessed from the ensemble of S matrices. In such a way one may obtain, e.g., the universal 
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conductance fluctuations lOI from the random matrix model fO . 



The study of mesoscopic devices points out that random matrix approach has its lim- 



itations, e.g., the fluctuations may become dependent on the length of the device ||15|| . It 
is still, an open question, what are the limits of the universality. We do not adress this 
problem here, rather we want to study, within the random model, the generic properties of 
S matrices dependent on some parameter. In this way we extend recent intensive studies of 
the statistical properties of bound systems dependent on the external parameters fI"6HI"9 . 



For parameter independent Hamiltonians statistical properties are generically universal 
once the mean level spacing, A is known. Parametric measures reveal a similar universality, 
in addition to unfolding the energy levels one has to unfold the parameter dependence |R|[17j • 
On the other hand, much less is known about the parametric behaviour of scattering systems. 
The two-point correlation function for the S matrix with respect to an external parameter 



was derived recently by Macedo [20]. The conductance fluctuations in the presence of the 
magnetic field induced time-reversal symmetry breaking have been studied by Pluhaf et al. 
PT[ . The semiclassical properties of the so called Wigner time delay have been discussed 



by Jalabert and Pichard [HJ. Obviously, even in the absence of the external parameter S 
matrix is energy dependent - the corresponding correlation function was obtained in 0. 

Instead of investigating the S"-matrix elements we focus on the properties of the cor- 
responding eigenphases (phase shifts). Some important dynamical features of the system 
(for instance the time delay inside the interaction region) can be easily expressed using 
the derivatives of S'-matrix eigenphases. Being directly accesible in experiments, statistical 
properties of phaseshifts deserve a detailed study. The nearest-neighbour spacing distribu- 



tion has been discussed in Ref. [p2p3| . As mentioned above we consider here parametric 
dependence of phaseshifts. 

The paper is organized as follows. The analytic results concerning the distribution of 
velocities (i.e., first derivatives of eigenphases with respect to energy or other external param- 
eter), also referred to as slopes, are discussed in Section II. Here we discuss also briefly the 
distribution of the second derivatives, i.e., the curvatures, of the eigenphases. These predic- 



tions are tested against numerical results obtained from simulations based on Eqs.( |l.l| - |1.2| ) 
(and their generalization allowing for the presence of some external parameter) in Section 
III. For systems dependent on the external parameter one may also construct directly cir- 
cular ensembles of scattering matrices. The corresponding results are presented in Section 
IV. 



II. PARAMETRIC HAMMILTONIAN APPROACH TO S MATRIX 



Consider the unitary matrix S defined by equations (|1.2|) and ( |1.3|) . To discuss properties 
of the eigenphases of the S matrix it is convenient to rewrite Eq. (|1.2p as 



with A given by 



S(E) = i±M (2 . 1} 



A = g 2 W + ^-^W. (2.2) 



Eigenphases s m of the unitary matrix S are related to eigenvalues a m of the M x M 
Hermitian matrix A by |§ 

s m = 2arctan(a m ) m — 1, . . . , M. (2.3) 

Since s m is a function of a m the statistical properties of s m and elm circ identical after 
unfolding. In the semiclassical limit (large M) it has been shown that A belongs to the 
same ensemble as H. This implies that the statistical properties (level spacing, number 
variance etc.) of s m are identical with those of the eigenvalues of H. Moreover, the relation 



( p.3| ) is useful when discussing the parametric dependence of s m . 

To see how it works consider the slopes of the eigenphases with respect to the energy E. 
They have a direct physical significance as time delays associated with the corresponding 
phaseshifts (see, e.g. |0|) while the average slope is just the celebrated Wigner time delay. 

The eigenequation for a m : 
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is equivalent to 



A\f m >= 9 2 W+^-^W\f m >= a m \f m > (2.4) 



H + — WW + ) \h m >= E\h m >, (2.5) 



with \f m > and \h m > related by \h m >= ^zjfW\f m >. Note that the N x N eigen- 
value problem Eq. (|2.5|) has, for fixed E, only M nontrivial solutions a m and corresponding 
eigenvectors \h m >. This is related to the fact that M — N eigenvalues of WW + vanish. 
Differentiating Eq.( [2.5| ) with respect to E for nontrivial a m one gets 

g 2 da m 



< h m \WW + \h m >=< h m \h m > . (2.6) 



a 2 n dE 

We shall assume from now on that \h m > are normalized to unity. 

Using the relation ( |2.3| ) between s m and a m , we express the inverse time delay as 

_i 1 1 + al 



in 



ds m / dE da m /dE 
Thus eliminating the derivative via Eq. (|2.6[) we obtain 



(2.7) 



u m = g 2 < h m \WW + \h m > +g 2 <h m \WW+\h m > _ m 



a 2 m 



But it follows from Eq. ( |2.4| ) and the relation between \f m > and \h m > that the second term 
on the right hand side above is proportional to the norm < f m \f m >■ Thus finally we get: 

u m = g 2 < h m \WW + \h m > + < f m \f m > /g 2 . (2.9) 

Eq.( [2.9| ) indicates that the distribution of inverse time delays is related to the norms ap- 
pearing on its rhs. Having in mind that < h m \WW + \h m > is a sum of M terms: 

< h m \WW + \h m >= | < h m \wi > | 2 + ... + | < h m \w M > | 2 , ( 2 -10) 

where Wk, k = 1..M, are vectors describing the coupling to the channel k, we can estimate 
the distribution of the matrix element < h m \WW + \h m > by a x 2 distribution with M degrees 
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of freedom in the GOE case and 2M degrees of freedom in the GUE case respectively. On 
the other hand < f m \f m > has a x 2 distribution with 1 (GOE case) or 2 (GUE) degrees of 
freedom. Assuming the two terms on the r.h.s of Eq.( |2.9| ) to be independent, we obtain the 
distribution P(u) by convolution of two xl distributions with different number of degrees 
of freedom and different means. The result can be expressed in terms of the confluent 
hypergeometric function iF\(a; c; x), sometimes called Kummer function |24|]. Changing the 
variable to r = 1/u we obtain the distribution of time delays 

P ( T ) _ expHfl 1 F 1 ^-^-tf- g -^ T ) 

{ } g 3(M-l) T ( ^+1) ^ r £^+l ' K ' 

where /3 — 1, 2 for orthogonal and unitary ensembles, respectively. 

Figure |2] shows the mean time delay (r) as a function of g obtained by integration of the 
above distribution for (3 = 2 and M = 2, 4, 6 and 16. In general mean time delay decreases 
with M. Moreover, the value of the coupling constant g m , for which the mean delay time is 
maximal, decreases with the number of channels. In both limiting cases g — > and g — > oo 
the mean time delay tends to zero, but the physical meaning of this fact is different. In the 
former case the coupling is so weak that the scattered wave is not affected by the bound 
system H. In the latter case very strong coupling causes the scattering to occur almost 
inst ant aneously. 

Formula ( |2.11|) simplifies in the special case g = 1. Then the inverse time delays are 
distributed according to x 2 distribution with f3(M + 1) degrees of freedom and the mean 
(u) = (3(M + l)/2. The time delay distribution, P(r) is then 

f 0-) = ^7B5J r ^ e ~' (2 ' 12) 

with the mean time delay 

(r) = 2/(/3(M + l)-2). (2.13) 

For large number of channels M the x 2 distribution can be approximated by a Gaussian 
with the variance equal to the mean (it). Thus for M large 



P(r) = -^L exp{-a 2 (2/r - 2/r ) 2 }, (2.14) 

where a 2 = l/(4/?(M + 1)) and r = 2//3(M + 1). 

Consider now the case when the Hamiltonian H = H(x) depends on some external 
parameter x and investigate the properties of s m (x) keeping the energy, E, fixed. 

Equation fl2.5|) , putting explicitly the x- dependence, takes the following form 

[H(x) + -£—WW + ) \h m {x) >= E\h m {x) >, (2.15) 

and may be viewed as an eigenvalue equation for eigenvalues E n (a, x) of the matrix H(x) + 
^WW + , defined for any real parameter a. Then Eq. fl2.15| ) is equivalent to a set of implicit 
equations 

E n {a,x) = E; n = l,2,..N (2.16) 

for an unknown a(x). Again due to the positivity of WW + the equation ( |2.16j ) has at most 
one solution a m (x) for each n with the total number of all possible solutions being equal to 
M. Moreover from Eq. fl2.16| ) follows 

^m(^) = dE m (a } x) c4 

dx dx g 2 < h m \WW+\h m >' { ' ' 

where E m denotes those eigenvalues for which ( |2.16|) has a nontrivial solution. Combining 
( p.l?l ) and (|273|) together with Eq. (|2^ ) yields simply 

d 9E m 

V m = -rS m = T m —— (2.18) 

dx ox 

providing the connection between the slope of the eigenphase with respect to external pa- 
rameter x, the corresponding time delay and the slope of the hermitian eigenvalue problem 
Eq. fl2.15p . The generic distribution of slopes for Hermitian random matrices is Gaussian-like 



|T7| , [r8"[ | . Thus the eigenphases slope distribution P(v ) may be obtained as a simple integral 
of the Gaussian and P(r) given by Eq. (|2.12|) [we consider the critical coupling, g = 1 case 



only]: 
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1 r°° ( ?; 2 \ d,T 

P(v) = ^= P(r)exp -— (2.19) 



2nJo v 7 \ 2r 2 / r 



under the assumption that appropriate unfolding of the parameter, x, has been made p7| , |18 
The integral in ( [2.19| ) yields |25j 



P ^ = 7/i^^ V -^y^M, (2.20) 



where T> p (z) is a shorthand notation for the product of the parabolic cylinder function, 
D p (z), |^5| with the exponential 

V p {z) = exp(z 2 /A)D p (z). (2.21) 

Using the series representation of T> p (z) for small z one easily verifies that the large \v\ tail 
of P(v) decays algebraically as P(v) oc l/|ti|^( M+1 )/ 2+1 . The asymptotic form of T> p (z) = 
z p [l + 0(z~ 2 )] |25j valid for large z yields regular behaviour of P{v) for \v\ small. 

Let us recall that the distribution of level slopes with respect to an external parameter, 
x, for a bounded generic system is Gaussian |T7| , |T8j . It is remarkable that the corresponding 
slope distibution for the eigenphases of the S matrix reveals the algebraic tails. As we shall 
demonstrate in the numerical example below, for large M the algebraic tail asymptotic form 
appears for very large \v\ only, while the center of the distribution resembles a Gaussian. 
One may verify, using the asymptotic form pf P(r) given by Eq.( p,14[) that the Gaussian 
distribution for P(v) is recovered in the M — > oo limit. 

A similar analysis can be carried out also for curvature distribution. In the limit of large 
M it can be shown in such a way that the curvature distribution tends to the generalized 
Cauchy one, as for the bounded systems [O . 



III. NUMERICAL RESULTS - HAMILTON! AN APPROACH 

We discuss the numerical results obtained for the slopes with respect to the energy, E, 
first. As discussed above they are directly related to possible time delays in the scatter- 
ing process and are, therefore, of particular interest. The time delay distributions, P(t), 
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obtained numerically by diagonalizing the S (|1.2|) for various number of channels M and 
g — 1 in the GUE case are presented in Fig.|3| and compared with the formula ( |2.12| ). The 



dimension of the internal Hamiltonian was taken iV = 100, energy E was set to zero and 
the data where obtained from 15000 generated 5* matrices. Fine agreement of the numerical 
results with the formula fl2.12|) does not depend on the energy and the number of eigenstates, 



provided N » M and E is not too large. Figure (f§) shows a comparison of the distribution 
P(t) obtained for g ^ 1 and the formula ( p.llj ). Again the agreement is quite satisfactory. 



Analogous results for the GOE case and g = 1 are plotted on the Fig.| (dashed line). 
The agreement with the prediction ( |2.12|) is not so nice as in the GUE case. Especially 
for small M a discrepancy is significant. Importantly, it may be removed when we use the 
formula ( [2.12| ) with one additional degree of freedom, i.e.: 



^ / n 1 M n 1 , , 

P ( T ) = ——— r --- 2 e --. 3.1 

- t M+2 ^ 



r 



2 

The corresponding average time delay is given by Eq.( [2.13| ) with (5 = 1 and M — > M + 1 
substitutions. We are not able to give a plausible explanation for this additive degree of 
freedom. The numerical results demonstrate clearly, however, that it is superior to Eq.( |2.12[ ) 
and leads to an excelent agreement with the numerical data - see Fig.[5| 

The formula (|3.1|) implies that the tail of the distribution for GOE behaves as r - M / 2 ~ 2 . 
In particular, for M = 3 we obtain for the time delay distribution P(r) ~ r _7//2 , in a full 
agreement with the experimental finding in ||26||. 



Similar numerical tests may be also performed for the velocity distribution in the case 
of a parametric dependence in order to test the formulae (|2.19| ),( [2.20| ). To this end we have 
to introduce a parametric dependence into the basic model of Eq.( |1.2| - |OD . One can, in 



principle, discuss several possible cases, the parameter x may affect either the bound system 
only [described by in Eq.( |1.3| )1 or the decay part, i.e. the bound-continuum coupling W, 
or both. In a generic case, arguably, both parts of the effective Hamiltonian will be affected. 
We consider the case when the internal dynamics only is x dependent, i.e., H = H(x). This 
situation corresponds directly to predictions obtained in the previous section. A similar 
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approach has been adopted in the treatment of time reversal symmetry breaking influence 



on the conductance fluctuations 21 



The parametric x-dependence is taken in a generic form: H — Hi cos(a;) + H2 sin(x) 
fT8| where Hi are drawn independently from the same ensemble of random matrices. The 
trigonometric form assures that the mean density of states remains the same for all x and the 



motion of bound levels as a function of x is stationary ||18|| . The slopes have been calculated 
by a finite difference in x from eigenphases of S. As before more than 15000 ^-matrices 
were used for the averaging. 

In the GUE case the agreement between the numerical data and the analytic prediction 
for the slope distribution, Eq. fl2.20|) is again remarkable - see Fig. || Note strongly non- 



Gaussian character of the obtained distribution. In the semi-logarithmic scale used in Fig. |6| 
a Gaussian would take an inverted parabola shape. The distributions show an algebraic 
tail P(v) oc ti~( M+2 ) for small M. For M large a center of the distribution resembles a 
parabola (i.e., a Gaussian distribution in the linear scale) followed by a straight line in the 
semilogarithmic plot indicating a regime of exponential behaviour, P(v ) oc exp(— ^y\v |). The 
numerical data are insufficient to detect a transition to a possible algebraic tail for \v\ large. 
For M even larger (not shown) the regime of Gaussian behaviour broadens, the exponential 
behaviour moves to the tails indicating the transiton to the semi classical limit. 

Consider now GOE internal dynamics. Similarly as for the time delays - the agreement 
with the numerical data is admirable if the distribution (|2.19|) is evaluated with the "ad- 



hoc" improved P(t) given by Eq.( |3.1|) - compare Fig. |^. As before the "proper" P(r) 
[Eq. (|2.12|) 1 with (3 = 1 used in arriving at Eq.( [2.20| ) leads to a theoretical prediction with clear 
disagreement with the numerical data for small M (not shown). The corrected expression 
is obtained from Eq. ( |2.20| ) by (3 = 1 and M — > M + 1 substitutions. The distributions 
obtained are similar in shape to those corresponding to GUE case, with the algebraic tail 
of the form P(v) oc t>~ M / 2 ~ 2 seen clearly for M small. For the largest number of channels 
plotted M = 16 we again observe the Gaussian center with exponential tails in similarity 
with the GUE distribution. 
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To complete the discussion of the parametric dependence in the Hamiltonian based ap- 
proach let us pass now to the curvature distributions. Numerical tests performed for the 
curvatures K = d 2 s/dE 2 are presented in FigM Note the quite good agreement of the 



numerical data with the generalized Cauchy distribution jrjj . The double-logarithmic scale 
used enhances the small and large \K\ behaviour. For large \K\ the agreement is excellent 
confirming the universality of the large curvature tail behaviour also for the present scat- 
tering system. On the other hand one observes a slight excess of small curvatures (and the 
corresponding lack of "medium" curvatures). It is a clear indication that M = 16 is not suf- 
ficient to realize fully the semiclassical limit. Using the analogy with the bound system level 
dynamics |18[ one may conclude that avoided crossings between eigenphases (as a function 



of the energy, E) are still partially isolated. 

The curvature distribution with respect to the external parameter is plotted in Fig.|9|. 
Comparing to the previous case, a better agreement with the generalized Cauchy distribution 
is observed indicating that the semiclassical regime is reached faster when the motion of 
eigenphases as a function of x is considered. A similar qualitative conclusion may be reached 
considering the distribution of time delays and of the slopes with respect to the parameter 
x. 

For M smaller one observes stronger deviations from the universal Cauchy distribution 
(not shown). The explicit dependence on the number of open channels in that case has 
obviously the same origin as the corresponding dependence in the case of slope (or time 
delay) distributions. 

IV. RANDOM PARAMETRIC 5-MATRIX 

As mentioned in the introduction one often employs random ensembles to directly model 
properties of S matrices [|K],|T^,|n| . One may envision a similar approach for the study of 



the parametric statistics. To this end one has to define the parametric dependence of S 
matrices directly, without reference to the underlying bounded dynamics. Obviously, there 



12 



is some ambiguity here, there are several possible choices. The ideal approach should be 
conceptually simple and, at the same time, lead to the same distributions as those obtained 
for the Hamiltonian based approach. After all in both cases the parameter independent 
matrices belong to generalized circular ensembles ||. We shall see below that the simplest 
possibilities agree with the Hamiltonian approach in the limit of large number of channels 
(semi classical limit) only. 

The first approach proposed takes S matrices in the form: 

S(x) = S exp(ixV) (4.1) 

where So is drawn from the appropriate circular ensemble (COE or CUE) of a given rank 
M while V is a Hermitian matrix (rank M), independent of Sq and drawn from GOE or 
GUE, respectively. Then the eigenphases at x = are the phase shifts of S . The slopes 
and curvatures of S are easily defined as the corresponding derivatives at x = 0. 

Such a form of a parametric dependence is "borrowed" from the typical form of a Floquet 
operator corresponding to quantum maps (e.g., for a famous kicked rotator or for the kicked 
top - for a detailed introduction to such problems see [27|). While in [^B] the form of the 



Floquet operator was determined by the dynamics of the system, here we assume a random 
Sq, as discussed above. 

Eq. (|4.1| ) implies a Gaussian distribution of slopes, since the slopes of eigenphases of S at 
x = are given by diagonal elements of matrix V represented in the eigenbasis of Sq. The 
latter are Gaussian distributed as we have taken Sq and V independently. The Gaussian 
character of the velocity distribution for the model ( f4.1|) is independent on the matrix size 
M and holds for both COE and CUE. Recalling the results obtained in the Hamiltonian 
model we see that the statistical predictions for eigenphase slope distribution coincide only 
in the M — > oo limit. 

The random S matrices belonging to a given ensemble of random unitary matrices (COE, 
CUE) may be generated (for arbitrary M) by drawing the generalized Euler angles from the 
appropriate probability distribution |29|| . Similarly one may construct the way to introduce 
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the parametric dependence by considering the infinitesimal changes of the generalized Euler 
angles [3D|. Our numerical results indicate that if one assumes that all M 2 — 1 angles 



are affected by the perturbation in a same way, the distribution of slopes is Gaussian for 
arbitrary M. 

Results received of 10 4 random unitary COE matrices with M = 40 are presented in 
Fig. [lO]. Gaussian character of the "velocity" distribution, obtained in this case, is not at all 
typical of any parametric dynamics defined for random unitary matrices. On contrary, one 
may construct many kinds of parametric dynamics in the space of unitary matrices leading 
to non-Gaussian distributions P(v). Exemplary data, represented in the figure by triangles, 
where obtained for a model in which only M diagonal Euler angles have beed varied. For 
small M the distributions obtained are closer to the results of the Hamiltonian model ( |1.2| ) 
than the pure Gaussian distribution. Similar results hold also for CUE case. 

Let us consider now the second derivatives of eigenphases with respect to x. For a Floquet 
operator corresponding to the kicked top model, the parametric curvature distribution has 
been studied already [^SJ and shown ]nj to obey the generalized Cauchy distribution with 



quite a good accuracy for all three universality classes. The numerical simulations performed 
by us indicate that for M large (semiclassical limit) the same property holds for the S matrix 
models described above. 

The simple propositions discussed above failed to reproduce the distribution of slopes 
with respect to the external parameter, x, obtained for small M within the Hamiltonian 
S'-matrix model. However this goal may be simply obtained by a slight modification of 
assumption concerning the matrix V in Eq.( |4.1| ). As discussed above the distribution of 
slopes of S matrix eigenphases is equivalent to the distribution of diagonal elements of V 
in the eigenbasis of So- Instead of choosing V to belong to the appropriate for the given 
symmetry Gaussian random ensemble one may chose V from the ensemble defined by the 
desired distribution of diagonal elements, Eq. ( |2.19[) . This defines the distribution of M 
elements of V leaving the remaining M(M — l)/2 elements undefined (in a chosen basis). 
This additional freedom may be utilized to make this ensemble less artificial than it seems 
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at first. While this aposteriori procedure may be hardly described as an elegant one, it 
gives an indication of the properties of the required ensemble of V matrices. In particular, 
following the discussion presented in Section II, it is clear that the construction of such 
an ensemble may involve scalar products of unnormalized vectors in M-dimensional space 
[compare Eq. ( |2.1(J| )]. Hopefully one may define the random ensemble for S matrices of the 



form (|4.1|) which is universal, i.e., reproduces not only the distribution of slopes but also 
of all other statistical measures involving the parameter x. This would allow to use the 
"direct" S matrix approach instead of the Hamiltonian one, also for the problems involving 
the external parametric dependence. 

V. CONCLUDING REMARKS 

In this work we have considered the statistical properties of derivatives of the S matrix 
eigenphases with respect to the energy, E, as well as with respect to some external param- 
eter x. Within the Hamiltonian approach we were able to derive, using simple heuristic 
arguments, the analytic distributions for the time delay valid for arbitrary number of open 
channels M and arbitrary value of the coupling constant g. An analytic distribution of 
slopes with respect to an external parameter has been also obtained. A comparison with 
the numerical simulations have shown good agreement of the theory with the numerical 
data in the CUE case. The discrepancies observed for the COE case may be removed by an 
"ad-hoc" modification of the proposed expression. 

When the main part of this work has been finished we have been informed that in the CUE 
case the same formulae for the time delay distribution have been obtained by Fyodorov and 
Sommers using supersymmetry calculus . We hope that the supersymmetric approach 



will be able in the future to verify the proposed distributions also for the time-reversal 
invariant (COE) model. 

We have verified numerically that the curvature distribution (i.e., the distribution of 
second derivatives of eigenphases with respect to energy) obeys the generalized Cauchy 
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distribution [18| in the semiclassical limit (M large). 

Construction of the parametric dependence directly for the random S matrices belonging 
to a given random ensemble has been shown to be ambiguous. The natural choices for the 
dynamics lead to Gaussian velocity distribution for arbitrary number of open channels. Then 
the agreement with the Hamiltonian based approach is obtained in M — > oo limit only. 
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FIGURES 

FIG. 1. Scheme of the scattering system. A system containing N bounded states and described 
by a Hamiltonian H is coupled via matrices W and W* to two waveguides with Mi and Mi 
open channels. The scattering in the system is thus characterized by a M x M matrix S, where 
M = Mi + M 2 . 

FIG. 2. Dependence of the mean time delay (r) on the coupling constant g for (5 = 2. Number 
of channels M labels each curve. 

FIG. 3. Time delay distributions, P(r), for GUE internal dynamics, g = 1 and for different 
number of open channels, M as indicated in the Figure. The numerical data (histograms) and 
the theoretical distributions (solid lines), Pi(t), Eq.( 2.12j ) are presented as a function of rj < r > 



to facilitate comparison of distribution shape for different M. According to Eq.(2.13) < r > is 
inversely proportional to M. 

FIG. 4. Time delay distributions, P(t), for GUE internal dynamics and g.ne.l: a) 
g = 0.5, M = 4; b) g = 2.0, M = 6. Solid line represents the theoretical prediction as given 
by Eg. (pip. 



FIG. 5. Time delay distribution, P(t), for time reversal invariant (GOE) internal dynamics, 
g = 1 and for different number of open channels, M as indicated in the Figure. The numerical 
data (histograms) are compared with both the analytical prediction, Eq. fl2.12D , (dashed lines) and 



the corrected prediction, Eq.(3T) (represented by a solid line) 



FIG. 6. Numerical data for distribution of slopes (with respect to external parameter) for GUE 
internal dynamic and g = 1 plotted in the semilogarithmic scale. The theoretical prediction, Eq. 
(|2.20j ) is shown as a full line and the number of channels is indicated in each graph. Note strongly 
non-Gaussian shape of the distribution. 

FIG. 7. Same as in the previous figure but for GOE internal dynamic. The theoretical predic- 
tion is obtained by the modification of Eq. Q2.20 ) discussed in the text. 
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FIG. 8. Distribution of second derivatives of eigenphases of S matrix with respect to energy, 
K = d 2 s/dE 2 obtained numerically for the case of M = 16 open channels, (a) for GUE internal 
dynamics, (b) for GOE dynamics. The solid (dashed) lines in both panels represent the generalized 
Cauchy distributions corresponding to the GOE (GUE) case respectively. 

FIG. 9. As in the previous figure for the curvature with respect to an external parameter, 
K = d 2 s/dx 2 . 

FIG. 10. Distribution of slopes of eigenphases with respect to the perturbation parameter x of 
random unitary matrices of size M = 40 typical of COE; a) parameter x controls variations of all 
M 2 — 1 Euler angles (A); b) only M diagonal angles are varied (o). 
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